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^ ■ Abstract 

^ ' In Phys. Rev. D 58, 014014 (1998) and 71, 094013 (2005), we determined non-perturbative 

, D*^ , Df, and fragmentation functions, both at leading and next-to-leading order in the 
MS factorization scheme, by fitting e^e~ data taken by the OPAL Collaboration at CERN LEPl. 
The starting points for the evolution in the factorization scale /i were taken to be /io = 2mQ, where 
. Q = c,b. For the reader's convenience, in this Addendum, we repeat this analysis for ^uq = mq, 

o 

' where the flavor thresholds of modern sets of parton density functions are located. 

o 

^ ; PACS numbers: 13.60.-r, 13.85.Ni, 13.87.Fh, 14.40.Lb 
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I. INTRODUCTION 



The OPAL Collaboration presented measurements of the fractional energy spectra of in- 
clusive D*^ D^, D^, and A+ production in Z-boson decays based on their entire 
LEPl data sample. Apart from the full cross sections, they also determined the contribu- 
tions arising from Z ^ bb decays. This enabled us, partly in collaboration with Binnewies, 
to determine lowest-order (LO) and next-to-leading-order (NLO) sets of fragmentation func- 
tions (FF's) for these charmed (Xc) hadrons 0, |^. We took the charm-quark FF to be of 
the form proposed by Peterson et al. and thus obtained new values of the e parameter, 
which are specific for our choice of factorization scheme. 

We worked in the QCD-improved parton model implemented in the pure modified 
minimal-subtraction (MS) renormalization and factorization scheme with nf = 5 mass- 
less quark flavors (zero- mass variable-fiavor-number scheme). This scheme is particularly 
appropriate if the characteristic energy scale of the considered production process, i.e., the 
center-of-mass energy ^/s in the case of e~^e~ annihilation and the transverse momentum 
Pt of the Xc hadron in other scattering processes, is large compared to the bottom-quark 
mass nib. Owing to the factorization theorem j^, the FF's defined in this scheme satisfy two 
desirable properties: (i) their scaling violations are ruled by the time-like Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) 7] evolution equations; and (ii) they are universal. Thus, 
this formalism is predictive and suitable for global data analyses. 

We verified that the values of the branching and average momentum fractions of the 
various c,b ^ Xc transitions evaluated at LO and NLO using our FF's in reasonable 

agreement with the corresponding results from OPAL P, IJ and other experiments j^. 

We tested the scaling violations of our D^, D^, D^, and A^ FF's by comparing the 
fractional energy spectra of these hadrons measured in non- resonant e~^e~ annihilation at 
y/s = 10.55 GeV 29 GeV Q, and 34.7 |ll| with our LO and NLO predictions to 
find reasonable agreement. Since events of Xc-hadron production from X^-hadron decay 
were excluded from the data samples at y/s = 10.55 GeV, we obtained a clean test of our 
charm-quark FF's. 

In Refs. jsl 0], the starting points /io for the DGLAP evolution in the factorization 
scale fi were taken to be fiQ = 2mQ, where Q = c,b. This choice is phenomenologically 
motivated by the observation that, in e"'"e~ annhilation, which has been providing the most 
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constraining input for the determinations of FF's, these values of /iq represent the very 
production thresholds of the respective flavors. Unfortunately, this choice is inconsistent with 
the convention underlying modern sets of parton density functions (PDF's) which prefer 
to place the flavor thresholds at /io = rnq. For the reader's convenience, in this Addendum 
to Refs. ^ I4I, we thus repeat the analysis of that papers for the choice /iq = rriq, so as to 
provide alternative LO and NLO sets of Xc FF's that can be conveniently utilized together 
with those PDF's. The FF's presented below were already used as input for a NLO study 
jisl ] of charmed-meson hadroproduction in pp collisions, which yielded agreement within 
errors with data collected by the CDF Collaboration in run II at the Fermilab Tevatron 
Q. we note passing that, in the case of pe.tu.bative.v .ndneed FF's. wh.ch ,s quite 
different from the case of non-perturbative FF's (involving substantial intrinsic components) 
under consideration here, the choice fio = mq is more natural, since, at NLO, it avoids finite 
matching conditions at the flavor thresholds jlS^ . 

II. RESULTS 

In the following, we concentrate on the most important results of Refs. ^, 0] that are 
affected by the shift in /iq. These include the fit parameters N , a, /?, and e defining the x 
distributions of the Q ^ Xc FF's Dq{x,^'^) at = /iq, 

Dt(x,^il) = Nx'H-xf, (2) 

the values per degree of freedom (x^/d.o.f.) achieved in the fits, and the branching 
fractions Bq{^) and average momentum fractions {x)q{^), 

Bgifi) = / dxDQ{x,fj.^), (3) 

J Xcut 
1 



{x)Q{f^) = p , . / dxxDQ{x,fi), (4) 

where Xcut = 0.1, at /i = 2/io and Mz- In the present analysis, we adopt the up-to-date 
input information from our 2005 paper Q . 

Our new results are presented in Tables HV lIVl Comparing Tables IIIII and IIVI with the 
corresponding tables in Refs. , we observe that the branching and average momentum 
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TABLE I: Fit parameters of the charm- and bottom-quark FF's in Eqs. ^ and ((21), respectively, for 
the various Xc hadrons at LO and NLO. The corresponding starting scales are fiQ = rric = 1.5 GeV 
and //Q = rrii, = 5 GeV, respectively. All other FF's are taken to be zero at = t^c- 



Xc Order Q N a j3 e 

LO c 0.694 0.101 

b 81.7 1.814.95 ••• 

NLO c 0.781 0.119 

b 100 1.85 5.48 • • • 

D+ LO c 0.282 0.104 

b 52.0 2.33 5.10 ••• 

NLO c 0.266 0.108 

b 60.8 2.30 5.58 • • • 

LO c 0.174 0.0554 

b 69.5 2.77 4.34 • • • 

NLO c 0.192 0.0665 

6 20.8 1.89 3.73 • • • 

Df LO c 0.0498 0.0322 

b 27. b 1.94 4.28 ••• 

NLO c 0.0381 0.0269 

b 27.5 1.88 4.48 • • • 

A+ LO c 0.00677 0.00418 

b 41.2 2.02 5.92 ••• 

NLO c 0.00783 0.00550 

b 34.9 1.88 6.08 • • • 



fractions are changed very little by the reduction in /xq- For a comparison of these observables 
with experimental data, we refer to Refs. 0, 0] . 

For lack of space, we refrain from presenting here any updated versions of figures included 
in Refs. j3[|4|; they would not exhibit any qualitatively new features. However, as already 
mentioned in Ref. [l3[, the reduction in /xq has an appreciable effect on the gluon FF's, which 
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TABLE II: x^/d.o.f. achieved in the LO and NLO fits to the OPAL ^, 3] data on the various 
hadrons. In each case, x^/d-o.f. is calculated for the Z ^ bh sample (6), the full sample (All), and 
the combination of both (Average). 
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NLO 


0.756 


0.560 


0.658 


D*+ 


LO 


1.19 


1.12 


1.16 




NLO 


1.07 


1.01 


1.04 


Dt 


LO 


0.246 


0.111 


0.178 




NLO 


0.290 


0.112 


0.201 




LO 


1.05 


0.117 


0.583 




NLO 


1.05 


0.112 


0.579 



TABLE III: Branching fractions (in %) of Q ^ Xc for Q = c,b and the various Xc hadrons 
evaluated according to Eq. Q in LO and NLO at the respective production thresholds = 2mQ 
and at the Z-boson resonance fi = Mz- 



X^ Order B^{2m^) B^{Mz) Bb{2mb) Bb{Mz) 





LO 


72.8 


67.6 


57.5 


52.7 




NLO 


71.6 


65.8 


54.3 


49.3 


D+ 


LO 


28.9 


26.8 


19.0 


17.7 




NLO 


26.4 


24.3 


18.5 


17.1 


D*+ 


LO 


29.0 


27.2 


24.3 


23.1 




NLO 


27.8 


25.9 


24.5 


22.8 


Dt 


LO 


12.3 


11.7 


23.1 


21.2 




NLO 


10.6 


10.0 


22.1 


20.2 




LO 


6.17 


6.06 


15.1 


13.7 




NLO 


6.12 


5.87 


14.3 


12.8 
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TABLE IV: Average momentum fractions of Q ^ X^. for Q = c,b and the various X^. hadrons 
evaluated according to Eq. @ in LO and NLO at the respective production thresholds = 2mQ 
and at the Z-boson resonance = Mz- 



Order {x)c{2m^) {x),{Mz) {x)b{2mh) {x)b{Mz) 





LO 


0.573 


0.442 


0.318 


0.285 




NLO 


0.550 


0.420 


0.304 


0.272 


D+ 


LO 


0.571 


0.441 


0.341 


0.302 




NLO 


0.557 


0.425 


0.324 


0.287 


D*+ 


LO 


0.617 


0.472 


0.393 


0.344 




NLO 


0.592 


0.448 


0.366 


0.322 


Dt 


LO 


0.654 


0.496 


0.348 


0.310 




NLO 


0.653 


0.487 


0.337 


0.299 


A+ 


LO 


0.765 


0.571 


0.302 


0.272 




NLO 


0.738 


0.544 


0.290 


0.261 



are only feebly constrained by e^e~ data. This effect is visualized for = D*^ in Fig. d 
where the /xq = mq to /io = 2mQ ratios of Dg{x,jj?) at /i = 5, 10, and 20 GeV are shown 
as functions of x. We observe that the reduction in /ig leads to a significant enhancement of 
the gluon FF, especially at low values of x. The results for X^ = D^, D^, D^, and are 
very similar and, therefore, not shown here. 



III. CONCLUSIONS 

In this Addendum to Refs. 3|, y], we repeated the fits of non-perturbative D^, D~^, D*~^, 
Dt, and At FF's, both at LO and NLO in the MS factorization scheme, to OPAL data from 
LEP: n a fo. the .educed choice ^ ^ c. .) of sta.t,„g point tot the DGLAP 

evolution in the factorization scale /x. These FF's are appropriate for use in connection with 
modern sets of PDF's |l2|, which are implemented with the same convention for the heavy- 
flavor thresholds. A FORTRAN routine that evaluates the values of these FF's as functions 
of the input variables x and fi may be obtained by electronic mail upon request from the 
authors. 
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This reduction in /iq is inconsequential for the theoretical interpretation of experimental 
e^e~ data because it is compensated by corresponding shifts in the fit parameters A^, a, /3, 
and e. However, the gluon FF's, which are only feebly constrained by e^e~ data, play a 
significant role in hadroproduction. In fact, detailed analysis [13] revealed that the increase 
in the gluon FF's due to the extension of the evolution length leads to a rise in cross section 
and thus improves the agreement with the CDF data of charmed-meson production in run II 



at the Tevatron 
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FIG. 1: /io = TUQ to /xo = 2mQ ratios of Dglxjfi"^) at /x = 5 (dashed), 10 (solid), and 20 GeV 
(dot-dashed) as functions of x for Xc = D*'^. 
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